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Abstract 



K*- ■ In this article, we consider an extension of the classical Curie- Weiss model in which 



the global and deterministic external magnetic field is replaced by local and ran- 
dom external fields which interact with each spin of the system. We prove a Large 



O ■ Deviations Principle for the so-called magnetization per spin Sn/n with respect to 



the associated Gibbs measure, where Sn/n is the scaled partial sum of spins. In 
particular, we obtain an explicit expression for the rate function, which enables an 
extensive study of the phase diagram in some examples. It is worth mentioning that 
the model considered in this article covers, in particularly, both the case of i. i. d. 
I random external fields (also known under the name of random field Curie- Weiss 

■ models) and the case of dependent random external fields generated by e. g. Markov 
chains or dynamical systems. 

■ Keywords: Disordered mean-field model. Large Deviations Principle, Random 

■ field Curie- Weiss model 
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1 Introduction 

Curie- Weiss models belong to the class of models in statistical mechanics for which 
one can explicitly explain important physical phenomena such as multiple phases, 
^ ■ metastable states and, particularly, how macroscopic observables fiuctuate around 

their m e an values when c lose to or at critical temperatures. Ellis and Newman (see 
EN78al . lEN78bl . IeNRSO^ ) computed the statistics of large spin-block variables in 



the class of classical Curie- Weiss models, which correspond to a finite collection of 
n (spin) random variables with an equal interaction of strength 1/n between each 
pair of spins, in the presence of a constant and deterministic external magnetic 
field. They derived a Law of Large Numbers and a Central Limit Theorem for Sn/n 
(where Sn is the total magnetization of the model with n spins), basically showing 
that Sn/y/n has a non-Gaussian limit at a critical temperature /3c, and a Gaussian 
limit otherwise. From a probabilistic point of view, a natural next step after show- 
ing a Law of Large Numbers and a Central Limit Theorem would be to study large 
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and moderate deviations propert ies oi Sn/ n. Wh ile a Large Deviations Principle 
(LDP) already goes back to Ellis [e1185I |. in jEL04l| Eichelsbacher and Lowe showed 
that Sn/n satisfies Moderate Deviations Principles (MDPs) with respect to the as- 
sociated Gibbs measure. 

In this paper, we are interested in the behaviom' of Sn/n for so-called Curie- Weiss 
models with random fields, which are derived from the Curie- Weiss model by replac- 
ing the constant external field by local and random external fields which interact 
with each spin of the system. Noteworthy exponents of these models are the random 
field Curie- Weiss models (abbreviated by RFCW models), where the local external 
fields are given by i. i. d. random variables, making RFCW models one of the easiest 
disordered mean-field models. These models have been studied i ntensively over the 
last decades, see e. g. jSPTTl . lAhaTi . ISWSSl . lAdMP9ll . lAdMPZ92l | and the references 



therein. It maybe worthwhile noting that many of the above authors assume the 
local external fields to be symmetrically Bernoulli distributed or to be bounded, 
while we do not need these restricting assumptions to consider RFCW models as 
representatives of Curie- Weiss models with ran dom fields . In the setting of RFCW 
models, Amaro de Matos and Perez showed in AdMP9l[ | a Central Limit Theorem 
for Sn/n, by using techniques similar to the ones employed by Ellis and Newman. In 
this fashion, they also proved that Sn/ \/n has a Gaussian or a non-Gaussian limit 
(depending on temperature), though the main difference with the classical setting 
is that these Gaussian and non-Gaussian limits depend on the realization of the 
local random external field. A natural next step, as before, would be to study an 
LDP and MDPs for Sn/n with respect to the now rand om Gib bs measure (disorder 
dependent). Almost sure MDPs have been proved in and an almost sure 

LDP can, in principle, be derived for RFC W models with bounded external fields 
using a result of Comets (Theorem V (ii) in Com89[ |). As hinted at, we can prove 
an LDP for RFCW models even if the external fields are unbounded. Yet, this is 
not the only reason why we want to prove an LDP in this article. The rate function 
obtained by using the result of Comets is extremly abstract and, hence, makes a 
study of the phase diagram impossible. We prove an LDP with an explicit rate 
function and make use of this advantage by studying the phase diagram in some 
examples. It should be mentioned that for certain choices of the random external 
field the analysis of the phase diagram for the corresponding RFCW model had been 
already obtained by using techniques related to the stud y of the dynamic s and the 
study of the metastates for the RFCW model ( see e.g . |FMPnnl . |bBIO^) and for 
more general mean-field models (see e.g. IK id . Kul97|). Moreover, our approach 



also yields an LDP if the external magnetic fields are not independent. 

The rest of the paper is organized as follows: In Section [2] we mathematically de- 
scribe the Curie- Weiss models with random fields. Our main result will be stated 
and proved in Section [3l Finally, Section S] is devoted to applications of our result 
and their use to study phase diagrams. 



2 Curie- Weiss models with random fields 

The Curie- Weiss model is a mean-field model of a ferromagnet. Due to its mean- field 
structure the spatial location of the spins is unimportant. The Hamiltonian of the 
Curie- Weiss model with external magnetic field /i G M can, therefore, be described 
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by 

^ n n 
i,j=l i=l 

(2.1) 

where Sn{cr) = Y17=i '^i total magnetization of the model. In the Curie- Weiss 

model with random fields the global external magnetic field h is now replaced by 
local and random external magnetic fields (hj,z G N), defined on some probability 
space {Q,J-,P), which satisfy P-almost surely 



1 " 

fn{x) ■■= - V'lncosh(x + /3hj) — > f{x) 

71 ^ ^ 



(2.2) 



n . 



as n — )■ oo for every x G Q and some differentiable function / : M — t- M. Its 
Hamiltonian is, thus, given by 

^ n n 

i,j=l i=l 

Correspondingly, the model can be associated with the following Gibbs measure on 
{ — at inverse temperature /3 > 

1 / /3 " " \ 

PnA^) = ^ exp ^ a.aj + /3 ^ h,c7, (2.3) 

where 

Zn,p = E I^E^^^^+z^E^^ 

o-e{-i,+i}" \ i,i=i «=i 
is the partition function of the model. 

We can see that, other than in the ordinary Curie- Weiss model, S„ is not an order 
parameter for the Curie- Weiss model with random fields, since in expression ()2.3p 
the measure is not completely determined by the value of 5„, now that the local 
fields also plays an important role. 

We want to end this section with three examples, which shall illustrate the diverse 
nature of Curie- Weiss models with random fields. 

Example 2.1 (Random field Curie- Weiss models). In this first example, the external 
fields are assumed to be i. i. d. random variables, that is, h is distributed according 
to the N-fold product measure of the marginal distribution of hi, i^. Assume that 
1/ has a finite absolute first moment, i. e. 




\h\ (iz/(/i) < oo. 

Since | lncosh(x)| < (|2.2p is, by the Law of Large Numbers, satisfied with 

f{x) = / lncosh(x -|- /3/i) (ii^(/i). 

JheR 

Example 2.2 (External fields generated by Markov chains). In this example, the 
external fields are still assumed to be identically distributed, but the independence 
assumption is dropped. To be more precise, let h be an irreducible Markov chain 
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on a countable state space S that has a stationary distribution vr. Moreover, let us 
assume that h is in equilibrium, that is, the initial state hi is distributed according 
to TT, and that vr has a finite absolut e first m oment. Then, with the help of the 
ergodic theorem (cf. Example 7.2.2 in Durld ]). ()2.2p is satisfied with 



f{x) = In cosh(x + ph) 7r(/i) . 

hes 

It should be noted that this function does not depend on the transition probabilities 
of the Markov chain, but only on the stationary distribution vr. 

Example 2.3 (External fields generated by dynamical systems). Let {0,,J^,P,(f)) be 
a dynamical system, that is, (17,7^, P) be the underlying probability space and 
(/) : Q — >• O be a measure preserving transformation. Assume that (p is ergodic 
and that X : — )• M is a J-"-measurable random vari able wi th finite absolute first 
moment. By the ergodic theorem (Theorem 7.2.1 in Durld l). ()2.2p is satisfied for 
(hi)ieN = (X o (p'-)ieNo with 



f{x) = [ lncosh{x + (3 h)dP^{h). 



3 Main result 

Before we state our result on an LDP for {Sn/n)n, let us recall the basic definition 
of an LDP. We will just consider the situation where elemen ts of i 3(M), the Borel 



(7-field on R, are of interest. For more general definitions see [DZ9i 



Definition 3.1. Let (P„)„gN be a sequence of probability measures. A sequence of 
real-valued random variables {Xn)neN is said to satisfy an LDP w. r. t. (Pn)nGN with 
rate function / : M — )■ [0, oo] if 

• L has compact level sets {x G M : L{x) < c} C M for any c G M, 

• (Upper bound) for every closed set C G B(^) it holds 

limsup-lnP„(X„ G C) < - inf /(x), (3.1) 

• (Lower bound) for every open set O G B{)E.) it holds 

liminf-lnP„(X„ G O) > - inf /(x). (3.2) 

n— >-oo n xGO 

Note that in this definition and the rest of this article the infimum of a function 
over an empty set is interpreted as oo. Also, what is called a rate function in this 
definition is sometimes referred to as a good rate function in literature. Having said 
this, we can state and prove our main result: 

Theorem 3.1. The magnetization {Sn/n)nem satisfies P-almost surely an LDP 
w.r.t. Pjj^ with rate function I{x) = J(x) — m.iy^^J{y), where J(x) = /*(x) — 
/3x^/2. Therein, f* denotes the Legendre-Fenchel transform of f, that is, 

f*{x) := sup{xy - /(y)}. 
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Proof. With probability one, we can choose a reahzation h = (/ij)jgf^ of h that 
satisfies ()2.2p . This reahzation satisfies ()2.2p not just for x £ but also for any 
x G M. Indeed, for x E M and every e > we can choose y G Q such that 

— y| < o and (3-3) 
\f{x)-f{y)\ < |, (3.4) 

since Q is dense in M and / is a continuous function. Note that for every n £ N fn 
is Lipschitz continuous with constant 1, since for xi,X2 G M, xi > X2, 

\fnixi) - fn{x2)\ " ^1 - X2I 



< 



. cosh(xi + /3/ii) 
cosh(x2 + f3hi) 



< 



^ ^ In ( ''vr^ : i -xi+x2 

i=l 

i V max (in f ^^!iii^l±M') , in ( ^0^X2 + ^h,) \ 1 ^ 
'^"^ I \cosh{x2 + (3hi) J Vcosh(xi +/3/ii)y J 

^ max In ^ ^3;j_|_3;2-(_y3/ij ^xi-X2-f}hi j \ pXi-X2+l3hi _|_ f>-xi-X2-l3hi ) | — 
1=1 



gXi+X2+phi _j_ QX\—X2—f}hi I ' \ g3;i— 3;2+/3/ii _j_ Q—x\—X2—f}h 

Therefore, (|3.3p implies 

l/n(x)-/„(y)| < |. (3.5) 

Since y G Q, we have \ fn{y) — f{y)\ < s/3 for n sufficiently large by assumption and, 
consequently, |/„(x)-/(x)| < |/„(x) - /„(y)| + |/„(2/) - /(y)| + |/(y) - /(x)| < e, 
where we have made use of ()3.4p and ()3.5p . Thus, we can assume the validity of 
(H^I) for any x G M. 

Next, let Qi for every i G N be the probability measure on { — 1, +1} induced by 



2 cosh(/3/ij 



and denote by Q the product measure 0°ZiQi- As a first step towards an LDP for 
Sn/n w. r.t. Pj^ p we prove an LDP w. r.t. Q , whic h will be done by means of the 
Gartner-Ellis Theorem (cf. Theorem 2.3.6 in j DZ98l |). To this end, we compute the 
logarithmic moment generating function A and observe that for x G M 

A(x) := lim llnEe"^^"/" 

n— >-oo n 

I Jl^ / gX+phi ^-x-/Shi\ 

= lim — > In ( I 

n-j-oo n \ 2cosh(pnj) / 

= lim (/„(x) - /„(0)) 

= /(^)-/(o), 

where E denotes the expectation w. r. t. Q. In doing so, we have used that Sn is 
a sum of independent random variables (under the product measure Q\) to derive 
the second line and that ()2.2p holds for any x G M to get the last line. Since / is 
differentiable, the same holds for A and the Gartner-Ellis Theorem yields that Sn/n 
satisfies an LDP w. r. t. Q with rate function R given by R{x) := f*{x) + /(O). Note 
that this implies f*{x) = 00 for |x| > 1 since {Sn/n)neN is uniformly bounded by 1. 
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Next, we see that o [Sn/n) ^ \s a. tilted version of Q o (Sn/n) ^, since with 

F -.R^R, X ^ {'^ ' II-' 
1/3/2, otherwise, 

and Xn = {—1, — 1 + 2/n, . . . , 1 — 2/n, 1} we have got for every S E B{R) 
Pll,po{SrJnr\S) = PnASn/n = m) 



I^meSnA-n e"-^*^™) H(T:S„{(T)=nm \\.i=l Qiii'^i}) 
I^meSnA-n e"-^(™') I]o-:5„(o-)=nm nr=l 

(m) 



The Tilted LDP (cf. Theorem III. 17 in dHOOl ]) then yields that Sn/n satisfies an 



LDP w. r.t. Pi's with rate function 

I(x) = R{x) - F{x) - mi{R{y) - F{y)} 
= nx)-F{x)-mf{ny)-F{y)} 
= r(x)-/3xV2-inf{r(y)-/3yV2}, 



where we have used f*{x) = oo for > 1 to obtain the last line, which, in turn, is 
the desired form of the rate function. □ 

Remark 3.1. Looking at the previous proof, we see that we carried over the problem 
to an LDP w. r. t. Q and proved this by means of the Gartner-Ellis Theorem. One 
could have, instead, thought about directly using the Gartner-Ellis Theorem to 
prove an LDP w. r. t. with rate function I. However, it is not only way more 
complicated to calculate the logarithmic moment generating function in that setting 
since Sn is not a sum of independent random variables under the measure P^^, but 
also does the logarithmic moment generating function in general not satisfy the 
assumptions of the Gartner-Ellis Theorem. This is evident, since the desired rate 
function / is in general not convex (cf. e. g. Example 14. 2p and the Gartner-Ellis 
Theorem can only be used to prove LDPs with convex rate functions. 

Remark 3.2. It is well-known from Large Deviations Theory (see e. g. Theorem II. 7. 2 
m [Eiisi ]) that the magnetization per spin Sn/n is asymptotically concentrated 
around the global minima of the rate function. Thus, it is important to know the 
precise structure of the rate function. The explicit representation of / given in 
Theorem 13.11 enables an accurate study of the global minima of the rate function 
/. Since / is the point wise lini it of convex functions {fn)n€N, f itself is convex and. 



thus, Theorem A.l in jCETOa ] yields that the global minima of / coincide with the 



global minima of G given by 



G{x) := Gp{x) := |x^-/(/3x). (3.6) 

Therefore, a study of G reveals the phase diagram of Curie- Weiss models with 
random fields and we identify the following phases of the system: 
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(i) Paramagnetic phase: G has a unique global minimum of type 1. 

(ii) Ferromagnetic phase: G has two global minima, both of type 1. 

(iii) First- order phase transition: G has several global minima, all of type 1. 

(iv) Second-order phase transition: G has a unique global minimum of type 2. 

(v) Tricritical point: G has a unique global minimum of type 3. 

In this process we said that a minimum m is of type A; G N if there is a positive real 
number A (called strenght) such that 

Gix) = Girn) + -T^ix -rnf^ + 0{{x -rnf^^^) as x ^ m. 
{2k)\ 

4 Examples 

In this section we want to apply Theorem 13.11 in some examples and study the 
corresponding phase diagram as elucidated in Remark 13. 2t 



4.1 Classical Curie- Weiss model 

Assume hj = /i for alH G N and some /i G M. In that case we have 

f{x) = lncosh(x + /3/i) 
and Theorem 13.11 yields that the magnetization satiesfies an LDP with rate function 

I{x) = sup {xy - f{y)} - ^x^ - ini \ sup {yz - f{z)} - ^y'^ 
Using the identity 

In cosh artanh X = — — ln(l — x^), 
which holds for all |x| < 1, one sees that 

sup {xy - f{y)} = Io{x) - /3hx, 

where Iq is given by 

/o(x) := 

Thus, we get 



ln(l + x) + ln(l - x) if|x|<l 
In 2 if|x| = l 

oo if Ixl > 1. 



/(x) = -^x'-phx + Ioix)-mf \-^y'-f3hy + Ioiy)}, (4.1 



which is the well-known representation of the rate function in the setting of the 
classical Curie Weiss model (see e.g. (4.17) in Ellsd ]). We omit a discussion of the 



rate function and the phase diagram as they are well-discussed (see e. g. (E1185t | 



7 



4.2 Curie- Weiss model with dichotomous external fields 

Assume that (hj)jgN is an i. i. d. sequence, where hi takes the values h and —h with 
probabihty 1/2 each, h £ M. Again, Theorem 13.11 yields that the magnetization 
satisfies an LDP with rate function / given by (cf. Figured]): 

1. If |x| < 1, then 

I{x) = ln2 - + |arsinh ^ (^6 + \/l + x^a^)^ + ^In ^ 

-- In (b + V^ + x^a^) - inf G(x) 
with a:=sinh(2/3/i),6: = cosh(2/3/i) and (cf. (ISlHjl ) 

G{x) = ^a;2-iln(cosh[/3(x + /i)]cosh[/3(x-/i)]). 



2. If Ixl = 1, then 



I(x) = In 2 - -x^ - inf G(x) 



3. If \x\ > 1, then 

I{x) = oo. 





Figure 1: I displays phase transitions as a function of /3 and h. 



By Remark 13.21 in order to obtain the phase diagram of the Curie- Weiss model 
with dichotoi nous extern al fields, it is enough to study G, which has been done in 
Chapter 5 of [AdMPZ92t . We note a rich family of phase transitions in /3 as well as 
in h: For h > ^ G has a unique global minimum at x = 0. In contrast to that the 
case h < ^ is more subtle. One finds a strictly increasing function / : [0, ^) — >■ M 
with /(O) = 1 and f{xn) — ?■ oo as x„ ^ such that 
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(i) P < f{h): 

G has an unique global minimum at x = 0. 

(ii) /? = f{h): 

G has an unique global minimum at x = if /i < | arcosh -y/3/2. Otherwise it 
has three global minima, one at x = and the others symmetric. 

(hi) /3 > f{h): 

G has two symmetric global minima. 

Consequently, we identify the following regions in the phase diagram: 

(i) If /i > ^ or /i < ^, /3 < /(/i), the system is in the paramagnetic phase. 

(ii) If /i < ^, /? > /(/i), the system is the ferromagnetic phase. 

(iii) The remaining region h < ^,(3 = f{h) is again delicate as that critical line 
consists of three different regimes. If /i < /ic = | arcosh y^3/2 we note a 
second-order phase transition. At h = he we note a tricritical point and for 
h > he = ^ arcosh y^3/2 we observe a first-order phase transition. 

The same analysis could have been done in the setting of non-symmetric dichoto- 
mous external fields, that is, hi takes the values h and —h with probabilities a and 
1 — a respectively, a G [0, 1]. In this case, we get 

G(x) = ^x^ — alncosh(/3(x + /i)) — (1 — a) lncosh(/3(x — /i)). 



obtaining the same phase diagram as the one presented by Theorem 4.1 in KLN07 |. 



4.3 Curie- Weiss model with uniformly distributed ex- 
ternal fields 

In this last example, we want to study the phase diagram of the Curie- Weiss model 
with uniformly distributed external fields, that is, let (hj)jgN be an i. i. d. sequence, 
where hi is uniformly distributed on the interval [—h, h] for some h > 0. Again, by 
Remark 13.21 we need to study G, which in this setting is given by 

= ln2 + ^x^-^^^°g(-^^P(-f' + '^')»+L(x) (4.2) 
2 4/3/1 ^ ^ ^ ' 




_§_(^2 I l2\ _ TT^ _ dilog(-cxp(-2/3(fe-|a;|))) i i ^ , 

~2h\-^ -I- J 24/3/1 4/3/1 u |X| ^ /t 

-;3|^| + dilog(-exp(-2/3(|x|-/.))) ii\x\>h. 

Therein, dilog denotes the dilogarithm that is defined by the power series 

diiog(z) = Y.^- 

n=l 

Very similar to what happened in the previous example, we now see (cf. Figure [2]) 
that there exists an increasing function / : [0, 1) — )■ M with /(O) = 1 and /(x^) oo 
as x„ 1 such that 

(i) /3 < f{h): 

G has an unique global minimum at x = of type 1 and, thus, the system is 
in the paramagnetic phase. 
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(ii) /3 > f{h): 

G has two symmetric global minima, both of type 1, and, thus, the system is 
in the ferromagnetic phase. 

(iii) /3 = f{h): 

On the critical line, the system shows a different behavior than the one dis- 
played in the previous example. In particular, we do not observe a first-order 
phase transition. 



Acknowledgements. The authors would like to thank Felix Lucka for an inter- 
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